


X 



Vortex-like Structures in the Skyrme - Einstein 

Chiral Model 

Yu. P. Rybakov, A. M. Tarabay, and I. G. Chugunov (Department of Theoretical 



'S ■ February 1, 2008 

^. 

Abstract 

> 

ps) ■ A regular method is suggested for constructing vortex-like solu- 

^«0 , tions with cylindrical symmetry in the Skyrme - Einstein chiral model. 

The method is based on the expansion of metric and field functions 
^^ ' in power series with respect to the two small parameters entering the 

^ ■ model. The length mass density of the vortex is estimated. 

^ . 1 Introduction 

r~| , The Skyrme chiral model [1] proved its efficiency in modelling the 

structure of baryons and light nuclei [2] via soliton-like configurations 
endowed with the topological charge (winding number), the latter 
being interpreted as the baryonic one. Recent years the interest to 

p^ \ general relativistic description of Skyrmions arose in view of possible 

astrophysical applications of topological solitons coupled to gravity 
[3-5]. The economy of the Skyrme - Einstein model in treating nu- 
clear and gravitational interactions of matter being its very attractive 
feature, one should however indicate the nonintegrability of the equa- 
tions of motion that forced the authors of the aforementioned papers 
to use numerical methods. In lieu of keeping up with this tradition 
we elaborate in the present work an analytic method of solving the 
equations for vortex configurations with the aim to approximate three- 
dimensional structures by closed vortices and to expand, in future, this 
method to the configurations with spherical or axial symmetry. 



2 Equations of the Skyrme - Einstein 
model: neutral vortex case 



The Lagrangian density of matter in the model in question can be 



constructed from the chiral current l,,. = U~^d^U, fi = 0,1,2,3, U & 



SU{2), and reads 
r - 



-^sp(y^) + ^sp([;^,z,]r,n), 



(1) 



where e, A are constant parameters. For the configurations with cyhn- 
drical symmetry the metric is chosen in the form 



ds' 



e^t^dt'' - e^^de - e^^d^^ - e^^dz^ 



(2) 



where the functions /U, a, (3, 7 depend on the radial coordinate ^ G 
(—00, +00). As a first step we consider the simplest case of the neutral 
vortex, when the chiral matrix U is taken in the form 



U = exp(ire(0), r 







—If 





(3) 



Putting (3) into (1) and taking into account of the Einsteinian 
gravity Lagrangian density Cg = R/2se, where R is the scalar curva- 
ture and ae is the Einstein gravitational constant, one finds the follow- 
ing expression for the action functional: 



A = 2tt fdt fdz fd^e''+^ 



ol^-^fR'J M ,,'R' 



ae 



1 






+ fi'^') - ^ ( e^-^@" + e"-/^ sin^ G 



(4) 



One infers from (4) that the theory admits the discrete symmetry 
/Li ^ 7. Therefore making the substitution ^ = 7 and imposing the 
harmonic coordinate condition [6] 



a = (3 + 2-/, 
one transforms the functional (4) as follows: 

A = 2tt fdt fdz fdS, 



(5) 



X (1 + 2e^X^e-^ sin^ 9) - ^e""^ sin^ G 



(6) 



where u = 4(/3 + 7), v = 4j3. 

The functional (6) admits the evident scale symmetry: 

^ —I- a^, z -^ az, u ^ u — 2lna, V ^ V, (7) 

that will be used later for simplifying the asymptotic behavior of so- 
lutions at space infinity. 

Now we introduce the dimensionless parameters in the problem by 
scale dilation: 

^ -^ ax, z -^ az, V ^ V + 2\na , u ^ u + \na , 

where a^ = 2A^. Then we come to the quasimechanical problem with 
X playing the role of " time" and with the Hamiltonian " action" (mea- 
sured in the units 2tt/)?) of the form 



S = - I dx 
2 



^ r /2 i2\ 
— (u — V ) 

V 



- e'2(l + 2e2e-^/2 sin2 6) - e""^ sin^ e] , (8) 

where v = 8ee/X'^. If one uses the natural units h = c = 1 and the 
experimental value of the parameter A = 2/Ft^, F-,^ = 186 MeV, then 
one gets v ~ 1.26 x 10"^^. The other dimensionless constant e appear- 
ing in the functional (8) takes the experimental value [7] e ~ 0.13. 
Thus there exist two small parameters in the problem: v ^ 10~^^ and 
e^ ~ 10~^, that allows one to utilize the perturbation scheme. 

First of all we write down the variational equations corresponding 
to (8): 



V 2 




(9) 


i^" = _ie«-- sin2 e - -e"^/2e'2 sin^ G, 
V 2 4 




(10) 


,2 
e"(l + e^e-''/^ sin2 6) - -g-^/^^'e' sin^ B = 






= -sin 26(6"-'' - e2e""/20/2)^ 




(11) 


and the "energy" integral 






i(t/'2 - ^;'2) + e«-- sin2 e = e'2(l + e^e"! sin^ 


6). 


(12) 



Using the scale symmetry (7) we impose the regularity conditions 
at space infinity: 

x^+oo, v^+oo, tt-u^O, e^O(e"''), (13) 



and on the axis [6]: 

X — > — OO, V — > — CX), ^ TT, 

I n — f |< CO, V exp I 1 -^ 4. (14) 

As follows from the equations (9) - (12), the regularity conditions (13) 
and (14) imply the asymptotic behavior of solutions at space infinity: 

V = u = 4:X, @ = yle~^, A = const, (15) 

and on the symmetry axis: 

V = Va + 4:Dx, U = Ua + 4:Dx, e = TT - Be^"", (16) 

with the integration constants Va, Ua, D, B satisfying the constraint 

Dew(^^^^^)=l- (17) 



3 Structure of vortex-like solutions 

Now we take into account of the mechanical analogy and rewrite the 
Eq. (12) in the form of Hamilton - Jacobi equation 

Sl = (l + e^e-^sin^e) x 

X [u{Sl - S^) + e"~^ sin^ e] . (18) 

Due to small value of the parameter i/ we can represent the solution 
to the Eq. (18) as the power series expansion 

S=-F + G + uH + ---. (19) 

Inserting (19) into (18) and equating the terms with the same powers 
in ly order by order, we get the equations for the unknown functions 
F,G,H,- ■ ■ . The equations for the function F read 

Fe = 0, (F„-F,)(F„ + F,)=0. (20) 

The asymptotic behavior (15) being equivalent to the relations (as 

X -^ oo): 

Fu = u' = v' = -F^, 

one concludes that the general solution to (20) turns out to be 

F = F{u-v). (21) 



The form (21) of the function F iniphes the following equation 
satisfied by the function G: 

G% = (l + e2e~isin2e)x 

X 2F'(u - v){Gu + G^) + e"-^ sin^ e] . (22) 

The variables in the Eq. (22) can be separated by the substitution 

F{u-v) = 8e("-^)/^ (23) 

oo 
n=0 

where the choice (23) is motivated by the nonuniform term in (22) 
and by the asymptotic behavior (15). Inserting (23) and (24) in (22) 
one obtains the recurrence relations satisfied by the functions Vn{Q)'- 

n 

J2 n'K_, + 4(n - 1) sin2 eK-i = 
fe=0 

= 6no sin O + Sni sin Q . (25) 

For n = one derives from (25) and (15) that 

Vb = 1 - cos e. (26) 

For n = 1 one gets, in view of (26), the equation 

2 sin eV( + AVi = sin^ 9. (27) 

Taking into account of the relations (15) and (16), one concludes 
that the solution Vi(G) to the Eq. (27) has the two branches Vi (B) 
with the domains 

D+ = {o<e<^},D- = {^<e< n}. 

The corresponding functions attain the forms: 



Ff = -1(1 + cos 6)2 (- — ^— + 3 - cose 
6 V 1 — cos d 



^1+ = ^(l-cose)sin2e. (28) 



For n > 2 one obtains the following integral representation for the 
solution to the Eq. (25): 

,_ n+cosoyj de /i-cose 
" ~~ u-cosej J sine Vi + cose 

a± 

X (^ E ^kVLk + 2(n - 1) sin^ QVn-ij , (29) 

where a+ = 0, a_ = vr.In particular, for n = 2 from (28) and (29) one 
deduces that 



i:il!m±2| +X:'4J(1-CO80)"), (30) 

n=l / 

where the following numerical coeficients are introduced: 
° = "135' 1 = 3' 2 = 9' ^3 = ^> 

Within the scope of the first two terms in the expansion (19) the 
following normal system of equations emerges: 

-n' = 5„ = -e 2 - + 2J«^ 2y„(e)], 

\ n=0 / 



1../ _ o _ lp(«-^)/2/^8 



zv 2 



+ ^e2"(n + l)e-"^/V„(G) , (31) 



n=0 



oo 
^ g{«-.)/2 ^ g2ng-nt,/2y^(0)_ 

n=0 

The system (31) should be complemented by the smooth matching 
condition at the point x = xq, where @{xo) = tt/2 : 



CO 



E e'"e-"^(-")/2 (v;+'(vr/2) - V^''{n/2)) = 0. (32) 

n=l 

In particular, within the confines of n = 2 approximation the Eq. (32) 
allows us to find, through the use of (28) and (30), the approximate 
value of the effective expansion parameter 

C = ie2e~"(^o)/2 ^ 1 ^ 0.232. (33) 

^ 3 161n2- 37/15 ^ ' 

Now we can use (33) to estimate the mass length density of the 
vortex in search defined as 

M=-^AS, (34) 

where AS* denotes the total variation of S{x) in the domains D : 
AS = 5+(oo) - S+{xo) + S-{xo) - S-{-oo) = 

= _2 - ie^^ £ 62"e-^(5+ - 5-)U=.o (35) 

n=l 

Using (31) one gets in the first approximation in z^ that 

u = v = 4x + 0{iy). (36) 

Inserting (36) into (35) and taking into account of (28) and (30), one 
finds up to C^ terms that 

A5«-2-4C-2C^ ("^-8^2^ (37) 

Substituting (37) into (34) one gets, in view of (33), the approximate 
value of the vortex mass 

An 
MRi -2-1.31. (38) 

A 

As for the structure of the function G(x) defining the radial distri- 
bution of matter inside the vortex, as well as those of the functions u{x) 
and v{x) describing the gravitational field, they can be derived from 
the system (31) through the use of the standard iteration technique. 

4 Discussion 

In conclusion it is worthwhile to underline that we considered here 
the case of the neutral vortex not endowed with the topological charge 



density. Therefore after closing the vortex piece we come to the config- 
uration with zero baryon number Q, this fact ensuing from the formula 
for the winding number 

Q = -^^"''' J ^'^ Mh h Ik). (39) 

However, this shortcoming can be cured by generalizing the substitu- 
tion (3) with the aim to generate the topological charge density in the 
vortex. To this end one can try the following ansatz: 

U = cos^ exp (icTsx) + ifJi sin 7/; exp (io"3(5), (40) 

where o"i, CJ3 are the Pauli matrices and ^, x, 5 stand for the chiral 
angles represented in the cylindrical coordinates p, (/3, z as follows: 

^ = i^{p), X = -kzr]{a- p), 6 = mip, (41) 

with k = const, m a Z. The radius a , entering the argument of the 
Heaviside step function rj, is chosen to satisfy the piece-wise smooth 
matching condition cosV'(a) = . In the simplest case we can impose 
the following boundary conditions: 

V^(O) = TT, ^{a) = ^, ^(00) = 0, (42) 

and also the closure one: 

/c/ = 27rn, neZ, z£ [-^,^], (43) 

where / stands for the length of the vortex piece. Inserting (40) into 
(39) yields, in virtue of (41), (42) and (43), the nontrivial baryon 
number 

Q = —^ (fx sin 2^ ([VV-Vx] V(5) = mn. (44) 

Thus, the suggested substitution meets all necessary requirements. 
The analysis of the Skyrme - Einstein vortex configurations endowed 
with the topological charge density will be exhibited in the subsequent 
paper. 



References 

[1] T. H. R. Skyrme, Nud. Phys. 31, 556 (1962). 

[2] R. M. Nikolayeva, V. A. Nikolayev, and O. G. Tkachov, Physics 
of Elementary Particles and Atomic Nuclei 23, 542 (1992). 



8 



[3] T. D. Lee, Phys. Rev. D35, 3637 (1987). 

[4] N. K. Glendenning, T. Kodanie, and E. R. Klinkhamer, Phys. 
Rev. D38, 3226 (1988). 

[5] M. Heusler, S. Droz, and N. Straumann, Phys. Rev. B285, 21 
(1992). 

[6] K. A. Bronnikov, Problems of the Theory of Gravitation and El- 
ementary Particles, No. 10, p. 37. Moscow: Atomizdat, (1979). 

[7] G. S. Adkins, Ch. R. Nappi, and E. Witten, Nucl. Phys. 228B, 
552 (1983) 



